Valley filter and valley valve in graphene 
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It is known that the lowest propagating mode in a narrow ballistic ribbon of graphene may lack 
the twofold valley degeneracy of higher modes. Depending on the crystallographic orientation of the 
ribbon axis, the lowest mode mixes both valleys or lies predominantly in a single valley (chosen by the 
direction of propagation). We show, using a tight-binding model calculation, that a nonequilibrium 
valley polarization can be realized in a sheet of graphene, upon injection of current through a 
ballistic point contact with zigzag edges. The polarity can be inverted by local application of a gate 
voltage to the point contact region. Two valley filters in series may function as an electrostatically 
controlled "valley valve" , representing a zero-magnetic-field counterpart to the familiar spin valve. 

PACS numbers: 73.23.Ad, 73.63.Rt, 85.30.-z, 85.75.Hh 



The potential of graphene for carbon electronics rests 
on the possibilities offered by its unusual band struc- 
ture to create devices that have no analogue in silicon- 
based electronics 111 Conduction and valence bands 
in graphene form conically shaped valleys, touching at 
a point called the Dirac point. There are two inequiva- 
lent Dirac points in the Brillouin zone, related by time- 
reversal symmetry. Electrons and holes in each valley 
are massless, with a large energy-independent velocity 
v = 10 6 m/s. Intervalley scattering is suppressed in pure 
samples 0, 0, • The independence and degeneracy of 
the valley degree of freedom suggests that it might be 
used to control an electronic device, in much the same 
way as the electron spin is used in spintronics [|| or quan- 
tum computing 

A key ingredient for "valleytronics" would be a con- 
trollable way of occupying a single valley in graphene, 
thereby producing a valley polarization. Different lines 
of approach are possible: One might use the fact that the 
lowest Landau level in a strong magnetic field is valley 
polarized 8, 9J. In zero magnetic field, the exchange in- 
teraction may lead to a spontaneous valley polarization 
[lfjj| . Here we propose an electrostatic method of valley 
polarization, that relies on the breaking of time-reversal 
symmetry by a transport current. 

Earlier work [II El El El G3 El 03 on S ra P hene 
ribbons (long and narrow ballistic strips) has shown that 
they may support a propagating mode arbitrarily close 
to the Dirac point, and that this mode lacks the val- 
ley degeneracy of modes that propagate at higher en- 
ergies. For armchair edges of the ribbon, this lowest 
propagating mode is constructed from states in both val- 
leys, but for zigzag edges only a single valley contributes 
[ill IT2I Il3l Il4 Il5| . In accord with time-reversal sym- 
mmetry, the mode switches from one valley to the other 
upon changing the direction of propagation. 

In this work we consider a quantum point contact 
(QPC), which is a short and narrow constriction with 
a quantized conductance G — n x 2e 2 /h 18]. (The fac- 
tor of two accounts for the spin degeneracy.) A current 



I is passed through the QPC by application of a volt- 
age difference V between wide regions on opposite sides 
of the constriction (see Fig. The orientation of the 
graphene lattice is such that the constriction has zigzag 
edges along the direction of current flow. We demon- 
strate by numerical simulation that on the first conduc- 
tance plateau the QPC produces a strong polarization of 
the valleys in the wide regions. This finding signifies that 
the valley polarization present inside the QPC is not de- 
stroyed by intervalley scattering at the transition from a 
narrow to a wide region. 

We show that the polarization of this valley filter can 
be inverted by locally raising the Dirac point in the region 
of the constriction, by means of a gate voltage, such that 
the Fermi level lies in the conduction bands in the wide 
regions and in the valence band inside the constriction. 
Two valley filters in series, one acting as a polarizer and 
the other as an analyzer, can block the current if they 
have the opposite polarity (see Fig. [3Jl. This extends 
earlier findings on the effect of a potential barrier in a 
graphene ribbon by Wakabayashi and Aoki Eil , demon- 
strating that a QPC can operate as a "valley valve" - 
a purely electronic counterpart of the magneto-electronic 
spin valve. 

Our calculations start from the tight-binding model of 
graphene, with Hamiltonian 
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The hopping matrix element 
and |j) are nearest neigbors on the honeycomb lattice, 
otherwise = 0. The electrostatic potential energy 
U{ = U(xi) varies only along the axis of the constriction. 
It equals Uoc in the wide regions and rises to Uo inside 
the constriction. We smooth the step wise increase of the 
potential over a length L s , according to the function 



if x < -L s /2, 
Qx.(ar) = { I + |sin(7ra;/L s ) if \x\ < L s /2, 

1 if x > L s /2. 
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FIG. 1: Schematic of the valley filter. Central panel: Hon- 
eycomb lattice of carbon atoms in a strip containing a con- 
striction with zigzag edges. Top panel: Dispersion relation 
in the wide and narrow regions. An electron in the first val- 
ley (modes n = 0, 1,2, . . .) is transmitted (filled circle), while 
an electron in the second valley (modes n = —1,-2,...) is 
reflected (open circle). Bottom panel: Variation of the elec- 
trostatic potential along the strip, for the two cases of an 
abrupt and smooth potential barrier (solid and dashed lines). 
The polarity of the valley filter switches when the potential 
height Uo in the constriction crosses the Fermi energy Ef- 




FIG. 2: Schematic of the valley valve (top) and corresponding 
potential profile (bottom). The current through the constric- 
tion is blocked if/iL = Ef — Ul and y,R = Ef — Ur have 
opposite sign. 
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FIG. 3: Dispersion relation of a graphene strip with zigzag 
edges. The spacing of the low-lying modes approaches A = 
^yfHitTa/W for W/a ^> 1. The zeroth and first modes have 
a larger spacing, approaching 3 A/2 for W/a ^> 1. Vertical 
lines mark the valley centers at k = 2n/3a and 47r/3a. 



The potential barrier U(x) — ?7oo + (Uo — £4o)[@l s (- t ) — 
Ql b (x — L)} is rectangular for L s = (solid line in Fig. 
n bottom panel), while it has a sinus shape for L s — L 
(dashed line). 

The dispersion relation of the honeycomb lattice in a 
strip with zigzag edges is shown schematically in Fig. 
H (top panel) and exactly in Fig. The wide regions 
support 2N + 1 propagating modes at the Fermi energy 
Ef, which form a basis for the transmission matrix t. 
Modes n = 1,2,... TV lie in the first valley [with lon- 
gitudinal wave vector ka £ (tt, 2tt)] while modes n = 
— 1, —2, ... — lie in the second valley [with ka £ (0, 7r)]. 
The zeroth mode n = lies in a single valley fixed by 



the direction of propagation. The conductance of the 
constriction is determined by the Landauer formula 

2e 2 N N 

G= —j— T n , T n = Knm| 2 - (3) 

n=-N m=—N 

The valley polarization of the transmitted current is 
quantified by 

p _ To + J2n=l ( T n - T -n) , .s 

C ~ rp ' ^ 

2^n=-N ±n 

where we consider the case (illustrated in Fig. that 
the zeroth mode lies in the first valley. The polarization 
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FIG. 4: Conductance (top panel) and valley polarization 
(bottom panel) for the valley filter of Fig. Q as a function of 
the electrochemical potential in the narrow region. The solid 
and dashed lines corresponds to abrupt (L s — 0) and smooth 
(L s = 8a) potential steps, respectively. 



P E [—1, 1], with P = 1 if the transmitted current lies 
fully in the first valley and P = — I if it lies fully in the 
second valley. 

We have calculated the transmission matrix numeri- 
cally by adapting to the honeycomb lattice the method 
developed by Ando for a square lattice [l^. Results are 
shown in Figs. 01 and We have fixed the width of the 
wide regions at Woo — 70\/3a (in units of the lattice 
spacing a). The electrochemical potential in the wide 
regions is set at Ep — Uqo = /ioo = r/3, corresponding 
to 2N + 1 = 29 propagating modes. The narrow region 
has width W = 20\/3a. We measure the electrochemical 
potential Ep — Uq = /io in the narrow region in units of 
the mode spacing A = ^^/STrra/W = 7rr/40 (indicated 
in Fig. El). 

The operation of the valley filter is demonstrated in 
Fig. 0] The top panel shows the conductance, while the 
bottom panel shows the valley polarization — both as a 
function of the electrochemical potential /io in the nar- 
row region. For positive /io the current flows entirely 
within the conduction band, and we obtain plateaus of 
quantized conductance at odd multiples of 2e 2 /h (as pre- 
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FIG. 5: Conductance for the valley valve of Fig.|21at fixed {il 
as a function of jiR. The solid and dashed lines correspond 
to abrupt (L a — 0) and smooth (L s — 8a) potential steps, 
respectively. 



dieted by Peres et al. [15|] )• Smoothing of the potential 
step improves the flatness of the plateaus (compare solid 
and dashed lines). The plateaus in the conductance at 
G = (2n + 1) x 2e 2 /h correspond to plateaus in the val- 
ley polarization at P = l/(2n + 1). On the lowest n = 
plateau, and for < /j,q < A, the polarization is more 
than 95%. 

For negative the current makes a transition from the 
conduction band in the wide regions to the valence band 
in the narrow region. This interband transition has pre- 
viously been studied in an unbounded system pfH l2l| . 
where it leads to selective transmission at normal inci- 
dence. In the QPC studied here we find that the inter- 
band transition destroys the conductance quantization 
— except on the first plateau, which remains quite flat 
in the entire interval — 3 A/ 2 < /io < 3A/2. The reso- 
nances at negati ve un a re due to quasi-bound states in 
the valence band [22ll23| . The polarity of the valley filter 
is inverted for negative /io, with some loss of quality (in 
particular for the smooth potential). 

The operation of the valley valve is demonstrated in 
Fig- El The current is blocked for — 3A/2 < /ir < with 
/ii on the first conductance plateau, so that the con- 
striction contains two valley filters of opposite polarity 
in series. The switching behavior of Fig. is similar to 
that obtained by Wakabayashi and Aoki in a simulation 
of a zigzag graphene ribbon containing a potential bar- 
rier |14j . Their geometry does not have well-separated 
valleys (since there is no wide multi-mode region), so an 
interpretation in terms of valley filtering is not as obvious 
as in the wide-narrow-wide geometry considered here. 

In conclusion, we have demonstrated by computer sim- 
ulation that a quantum point contact with zizag edges is 
capable of producing a highly nonequilibrium population 
of the valleys in a carbon monolayer. The valley polar- 
ization is better than 95% if the Fermi level lies above 
the Dirac point both inside and outside of the constric- 
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tion. The polarity of the valley filter is inverted when 
the Fermi level in the wide and narrow regions straddles 
the Dirac point, permitting the construction of a switch 
(valley valve) consisting of two filters in series. We an- 
ticipate that the experimental realization of the device 
proposed in this work will make it possible to exploit the 
valley degree of freedom, in addition to spin and charge 
degrees of freedom, as a carrier of information in carbon 
electronics. 
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